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Abstract 

We derive for deep-inelastic neutrino(v)-proton(P) scattering in the combination vP — vP the per- 
turbative QCD corrections to three loops for the charged current structure functions F2, El and i% 
In leading twist approximation we calculate the first five odd-integer Mellin moments in the case 
of F2 and Ft and the first five even-integer moments in the case of F3. As a new result we obtain the 
coefficient functions to o(a? s ) while the corresponding anomalous dimensions agree with known 
results in the literature. 



1 Introduction 



Predictions for structure functions in deep-inelastic scattering (DIS) including perturbative cor- 
rections in Quantum Chromodynamics (QCD) have recently been advanced to an unprecedented 
level of precision over a wide kinematical region of Bjorken x and Q 2 = —q 2 , with q being the 
momentum of the exchanged gauge boson. The knowledge of the complete three-loop splitting 
functions for the scale evolution of unpolarized parton distributions of hadrons [1,2] together 
with the second-order coefficient functions [3-7] has completed the next-to-next-to-leading order 
(NNLO) approximation of massless perturbative QCD for the DIS structure functions F\, F2 and 
F3 . In addition for electromagnetic (photon-exchange) DIS the three-loop coefficient functions for 
both F2 and Fi = F2 — 2xF\ have become available [8,9], the latter being actually required to com- 
plete the NNLO predictions, since the leading contribution to the coefficient functions is of first 
order in the strong coupling constant a s . 

In the present article, we extend the program of calculating higher order perturbative QCD 
corrections to the structure functions of charged current DIS. Our studies are motivated by the 
increasingly accurate measurements of neutral and charged current cross sections at HERA with a 
polarised beam of electrons and positrons [10-12]. At the same time we are also able to quantita- 
tively improve predictions for physics at the front-end of a neutrino-factory, see e.g. Ref. [13]. To 
be specific, we consider neutrino-proton scattering in the combination vP — \P, which corresponds 
to charged lepton-proton DIS as far as QCD corrections are concerned. Following Refs. [14-18] 
we compute the perturbative QCD predictions to three-loop accuracy for a number of fixed Mellin 
moments of the structure functions F2, Fi and F3. 

Within the framework of the operator product expansion (OPE), and working in Mellin space, 
pvP-vP an( j pvP-vP ^ f unct j ons Q f Mellin moments only, while only even moments con- 
tribute to F^ ~ . This distinction between odd and even Mellin moments is opposite to the case 
of the neutral current structure functions and also to the case of charged current structure func- 
tions for neutrino-proton scattering in the combination \P + \P. In the latter case, the three-loop 
results for F 2 V P+vP and fV p + vP C an be directly checked in electromagnetic DIS and taken over 
from Refs. [8,9]. Also F^ p+vP is known to three-loop accuracy [19] with parametrizations for the 
respective coefficient functions given in Ref. [20] . 

Having available a limited number of fixed Mellin moments for F2, Fi and F3 for both com- 
binations of neutrino-proton scattering, i.e. \P±vP is a prerequisite for a subsequent complete 
calculation of the respective quantity to three loops. With the methods of Refs. [1,2, 8,9] at hand 
we have all ingredients for a future computation of the "all-n" results in Mellin-n space, or equiv- 
alently the complete expression in Bjorken-x space after an inverse Mellin transform. However, 
applying the present results we can already comment on a number of phenomenological issues, 
which we do in a companion paper [21]. 

The outline of this article is as follows. In Section [2] we briefly recall our formalism, which 
is based on the optical theorem, the forward Compton amplitude and the OPE. Specifically we 
emphasize the symmetry properties of the Compton amplitude for neutral and charged current pro- 
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cesses and show how these select either odd or even Mellin moments for the structure functions 
F2, Fi and F3 depending on the process under consideration, i.e. \P±vP. In Section [3] we recall 
details of the renormalization and give all relevant details of the calculation in Section HI Sec- 
tion [5] contains our results for the Mellin moments of F^ v , F^ _v and F^ ~ v in numerical 
form. Finally, we conclude in Section [6l The lengthy full expressions for the new moments of 
the coefficient functions are deferred to Appendix [A] and some details on the OPE are given in 
Appendix [Bj 



2 General formalism 

The subject of our calculation is unpolarized inclusive deep-inelastic lepton-nucleon scattering, 

l(k) + nucl(p) l'(k')+X, (2.1) 

where l(k), l'(k') are leptons of momenta k and k', nucl(p) denotes a nucleon of momenta p 
and X stands for all hadronic states allowed by quantum number conservation. In this article we 
are concentrating on charged current neutrino(v)-proton(/ 5 ) scattering, i.e. vP, \P via W ± boson 
exchange. As it is well known, the differential cross section for reaction (12.11) can be written as a 
product of leptonic and hadronic Wu V tensors 

da oc ZTWuv . (2.2) 

The leptonic tensor IF 1 for electroweak or pure electromagnetic gauge boson exchange is detailed 
in the literature, see e.g. Ref. [22] and will not be considered here. The hadronic tensor in Eq. (12.21) 
is given by 

W^{p,q) = d\^ z {nuc\,p\jl{z)MO)\nnc\,p) (2.3) 

= ^F L (x, Q 2 ) + d MV j-F 2 (x, Q 2 ) + ie^ a p|-^ F 3 (x, Q 2 ) , 

where is either an electromagnetic or a weak current and |nucl,p) is the unpolarized hadronic 
state with momentum p. The boson transfers momentum q, Q 2 = —q 2 > 0, and the Bjorken scaling 
variable is defined as x = Q 2 / (2p-q) with < x < 1 . The tensors and d^ v are given by 

e^v = g/iv-— 2~, (2-4) 

4x 2 2x 
dps, = -g/jv-PuPv— -(Pnqv + Pvqfi)^, (2.5) 

and £^ V(X p is the totally antisymmetric tensor. The hadron structure functions F ; , i = L, 1,2,3 are 
the main subject of our investigations in the present paper, with F\ being related to Fi and F2 by 
the Callan-Gross relation, 

F L (x, Q 2 ) = F 2 (x, Q 2 ) - 2xF x (x, Q 2 ) . (2.6) 
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The structure function F3 describes parity-violating effects that arise from vector and axial-vector 
interference and vanishes for pure electromagnetic interactions. 

We are interested in the Mellin moments of the structure functions, defined as 

1 



Fi(n,Q 2 ) = Jdxx'^F^Q 2 ), i = 2,L; (2.7) 


1 

F3(n,Q 2 ) = f dx^- l F 3 (x,Q 2 ). 



(2.8) 



The optical theorem relates the hadronic tensor in Eq. (12.31) to the imaginary part of the forward 
scattering amplitude of boson-nucleon scattering, 7^ v , 

W M v(p,q) = ImT^v (/?,#). (2.9) 

The forward Compton amplitude T pV has a time-ordered product of two local currents, to which 
standard perturbation theory applies, 

T, v (p,q) = i|j 4 ze 19Z (nucl,p|r(7;(z)7 v (0))|nucl,p). (2.10) 

In the Bjorken limit, Q 2 — > °°, x fixed, the integral in Eq. (12.101) is dominated by the integration 
region near the light-cone z 2 ~ 0. In this region the phase in the exponent in Eq. (12.101) becomes 
stationary for the external momentum q being deep in the Euclidean region. Thus, we can use the 
OPE for a formal expansion of the current product in Eq. (12.101) around z 2 ~ into a series of local 
composite operators of leading twist (see e.g. Ref. [23] for details). In terms of local operators for a 
time ordered product of the two electromagnetic or weak hadronic currents the OPE for Eq. (12.101) 
can be written in the following form 

i J J 4 ze^r(yt(z)7 v (0)) = 2£(^)"[(^ V -^)^^ 2 C L ,;(^^^) (2.H) 



Q 2 

-(g m gvn 2 q 2 - g m qvqn2 -gv^q^q^ +g f jvq^q^c 2 j(n, %,a 5 ) 



n,j 

e 2 

/ O 2 \ 1 

+ie Wl p/V« c 3,j (n, ^,a s ^q^...q^O j ^--^{^) + higher twists, 

where j = a, q,g and all quantities are assumed to be renormalized, /j being the renormalization 
scale. Higher twist contributions are omitted in Eq. (12.1 II) as they are less singular near the light- 
cone z 2 ~ and suppressed by powers of l/Q 2 . Therefore, the sum over n in Eq. (12.1 II) extends 
to infinity and runs only over the standard set of the spin-n twist-2 irreducible symmetrical and 
traceless operators. In a general case three kind of operators contribute (these correspond to the 
index j in Eq. (12.1 II) ): the flavor non-singlet quark operators O a , the flavor singlet quark operator 
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O q and the flavor singlet gluon operator O g . These are defined by, 

a,{n u -,»n} = vj/^yiw Z^'V, a= 1,2,..., (4-1), (2.12) 
ObiMu-'M = vj7y{wzy 2 ---Zy"V 5 (2.13) 

0g>{w,-,M _ pV{mjy*2 . . ,£)H„-lpMn}V (2.14) 

Here, \\f defines the quark operator and F^ y the gluon operator. The generators of the flavor group 
SU (n f) are denoted by X a , and the covariant derivative by D M . It is understood that the symmetrical 
and traceless part is taken with respect to the indices in curly brackets. 

The spin averaged operator matrix elements (OMEs) in Eqs. (I2.12I) - (I2.14I) sandwiched between 
some hadronic state are given by 

(nucl^O^ 1 ' -'^Inucl,/?) = P {mi ■■■P fJ " } A{ ucl (n,fj 2 ) , (2.15) 

where hadron mass effects have been neglected. The OMEs themselves as given in Eq. (12.151) are 
not calculable in perturbative QCD, but they can be related to the quark and anti-quark distributions 
of a given flavor and to the gluon distribution in the hadron. 

The scale evolution of the OMEs governed by anomalous dimensions is accessible to per- 
turbative predictions as well as the coefficient functions Q ; multiplying the OMEs according to 
Eq. (12.1 II ). Both the anomalous dimensions and the coefficient functions are calculable order by 
order in perturbative QCD in an expansion in the strong coupling constant a s . In order to do so, 
we replace the nucleon state |nucl,p) in Eqs. (12.101) . (12.1 II) by partonic states. In complete analogy 
to Eq. (12.101) we define the forward Compton amplitude fa at parton level and the corresponding 
partonic OMEs, 

(parton,p|0^ 1 '-"^ } |parton,p) = p^.-.p^A^^ (n./u 2 ) . (2.16) 

As the OPE in Eq. (12.111) represents an operator relation, we derive the following parton level 
expression 

fa = i /j 4 ze i9z (parton,/7|r^(z)7 v (0)) |parton,/>) (2.17) 

= 2 £ a>" C LJ in, — , a s J + C 2 j [n, — , a s j 
p a qfi / Q 2 \-i 

^parton { n i A' 2 ) + higher twists , 



which is an expansion in terms of the variable (0 = (2p ■ q)/Q 2 = 1 fx for unphysical (0 — > (x — > °°). 
The coefficients Qj with i = 2, 3,L are of course the same as the previous ones appearing in 
Eqs. (12.1 II) and the scale evolution of the OMEs in Eq. (12.161) is controllable in perturbation theory. 

Let us in the following recall a few aspects of flavor (isospin) symmetry of the DIS structure 
functions which are relevant to neutrino-nucleon scattering. The composite operators Eqs. (12.121) - 
(12.141) are either singlet or non-singlet operators referring to the representation of the SU (n/) flavor 
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group. In particular the non- singlet operator O a ^ 1 "''^ in Eq. (12.121) contains the generators X a 
of the flavor SU («/). It is well known, that for the separation of the singlet and non-singlet contri- 
butions to structure functions, Wilson coefficients, etc., one considers the sum and the difference 
of matrix elements for a proton P and a neutron N, e.g. 

F eP±eN = F eP ±F eN ? p vP±vN = f vP ± p vN ^ / = 2,3,L. (2.18) 

The combination P + N singles out contributions to the singlet (isoscalar) operators and P — N the 
corresponding ones to the non-singlet (isovector) operators, which can be seen readily as follows. 
To that end, let us specialize for simplicity to the case of two flavors (rif = 2) only, i.e. to a SU (2)- 
isospin symmetry, the generalization to an arbitrary number rif of flavors being straightforward. 
Then, in the SU (2) example, the twist-two term of the OPE consists of an isoscalar (Go) and 
isovector (0 a ) part, i.e., 

= o l + e a X a , a= 1,2,3, (2.19) 

where 1 is unit matrix and A, a = c a /2, o a are the usual Pauli matrices in fundamental representa- 
tion. Sandwiching Eq. (12.191) between the proton \P) and neutron \N) states, one gets directly 

(p\&\p) + (n\&\n) = e (p\p)+e a (p\i a \p)+%(N\N) + e a (p\i a \p) (2.20) 

= 90 + ^93 + 90-^93 = 290, 

(P\©\P)-(N\®\N) = 9 + ^93-(9o-^9 3 )= 93. (2.21) 

Here one uses the fact that proton and neutron are eigenvectors of the X3 isospin operator with 
eigenvalues +1/2 and —1/2, respectively. Hence the combinations of the OMEs such as A J p±N (n) = 
Ap(n) +A J N (n) correspond to the isoscalar part (singlet contribution) and isovector part (non-singlet 
contribution), respectively. As an upshot one can conclude that the OPE for the P — N combination 
receives contributions from non-singlet quark operator O a Eq. (12.121) , i.e. j = a, in the r.h.s of 
Eq. ( 12.1 II) . On the other hand, for the P + N combination both singlet quark operator O q Eq. (12.131 ) 
and singlet gluon operator O g Eq. (12.141) contribute in the OPE, i.e. the sum on the r.h.s Eq. (12.1 II ) 
runs over j — q, g. 

Since in the present article, we are considering charged current DIS in the combination \P — vP 
we have due to isospin symmetry 

^ F vP_ F vP = F vP_ F vN = F W_ F vN (222) 

Thus, we are entirely restricted to non-singlet quark operators for the structure functions i 7 ^ - ^, 

Next, we would like to address the symmetry properties of the partonic forward Compton 
amplitude (12.171) t^ v and explain how these translate into selection rules for either even or odd 
Mellin ^-moments of the different DIS structure functions. To that end, let us inspect the Feynman 
diagrams for at the leading order with initial quarks in Fig.Q] There the right diagram is simply 
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Figure 1 : Leading order diagrams contributing to the forward Compton amplitude in deep-inelastic 
boson(V)-quark scattering. 

the crossed diagram of the left one. In Fig. Q] we denote the gauge bosons by and V v . For the 
latter there are various choices as they can either be photons y or weak gauge bosons Z° and W ± . 
The matrix element of the left diagram is proportional to 

t^^T v p(p + q)T^ (2.23) 

where and r v denote the vertices of vector boson-fermion coupling, while p(p + q) is the quark 
propagator of momentum p + q, p(p + q) = — \/{j> +4)- For the right diagram one has 

t^ocT^pip-q)^. (2.24) 

Under the simultaneous transformation /j <-> v and q — > — q the matrix element of the crossed 
diagram is equal to uncrossed one, provided both vertices 1^, and T v have the same structure. 

Let us detail this situation for the case of the neutral current DIS first. The external bosons 
and Vy in Fig.[T]being photons y or Z° -bosons couple to the vertices and r v . The latter are either 
proportional to e q y^ and e q y v with the fractional quark charge e q (y-boson) or to (v/fy + 0/7^75) 
and (v/Yv + <3/YvY5) with the (flavor-depended) vector and axial- vector current coupling constants 
Vf and af (Z° -boson). In the case of y— Z° -interference one has to consider both, y and Z , in the 
initial state in Fig. \T\ with a different gauge boson in the final state. In the end the effective number 
of diagrams for the interference contributions will be doubled. For all neutral current DIS cases 
the quark flavor, of course, remains conserved. 

At this point, we can relate the action of simultaneously transforming /j <-> v and q — > —q 
in all Feynman diagrams contributing to t^ v to the parameters of the OPE in Eq. (12.171) . namely 
the coefficient functions C2, C3 and Cl- It is clear that the full matrix element for (l.h.s. of 
Eq. (12.171) ) remains unchanged, since the transformation /j <-> v and q — > — q maps the crossed and 
uncrossed diagrams onto each other, even in the case of y — Z°-interference due to the doubled 
number of diagrams. On the r.h.s of the OPE in Eq. (12.171) the tensors and d^ v remain invariant 
under /j <-> v, while the antisymmetric tensor e^, va ^ picks up a sign (-1). The coefficients C2, C3 
and Cl as well as the OMEs A part on, being Lorentz scalars, are at most functions of Q 2 = —q . 
Therefore they are invariant as well. Finally (0 will be transformed to —CO (recall its definition 
00 = (lp-q)/ Q 2 ). Thus, in the series expansion in spin n in Eq. (12.171) the coefficient functions C2 
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and Ci are weighted by a factor ( — 1)", whereas C3 is multiplied by ( — 1)" . In other words, the 
sum in Eq. (12.171) runs for C% and Cl only over even Mellin moments n and only over odd for C3. 
The coefficients for other n have to vanish in Eq. (12.171) . The same choice of n holds for the Mellin 
moments of the structure functions F2 and Fi (even n) Eq. (12.71) and F3 (odd n) Eq. (12.81) of neutral 
current DIS because of relations Eq. (12 .3 0112.311) which will be discussed later. 




Figure 2: Representative three-loop diagrams for the various flavor classes in charged current 
neutrino-proton DIS (see text). 

So far our discussion has been based on leading order Feynman diagrams but the previous 
arguments carry over to higher orders as well. Up to three-loop accuracy all Feynman diagrams 
fall in one of the so-called flavor classes fh, flo2, fhi or //on displayed in Fig.[2l The class fh 
corresponds to all diagrams with both gauge bosons and V v attached to the open fermion line 
of the initial (final) state quark. Class /Z02 collects the diagrams with both gauge bosons attached 
to an internal closed fermion loop, while fl\\ contains the diagrams with one boson attached to 
the closed loop and the other to the open line of the external quark. Finally the class //on denotes 
diagrams with both bosons attached to different closed quark loops. Depending on the process 
under consideration some flavor classes vanish identically. It is easy to see that the neutral current 
DIS assignments for C2 and Ci (even Mellin moments) and C3 (odd Mellin moments) from above 
persist and the same holds true for the structure functions. 



Figure 3: Leading order diagrams contributing to the forward Compton amplitude of charged 
current vP and vN scattering. The right diagram represents the crossed of the left one but with an 
incoming quark of different flavor. 

Let us next turn to the case of charged current DIS. We have the structure functions F2, F3 and 
Fl for both, an isoscalar and an isovector target, i.e. vP±\/N, which we have to distinguish (see 
also Eq. (12.221) ). On the partonic level this implies that we sum the contributions of u and d quarks 
in the singlet case and take their difference in the non-singlet case. For simplicity, we restrict 
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ourselves here again to SU (2)-isospin symmetry with flavors u and d only. The generalization to 
s, c and more flavors should be clear. In charged current vP or vN DIS we are considering initial 
and final gauge bosons V M = and V v = Wy + (cf. Fig. [T]), the coupling of J-quarks to W~ being 
excluded by electroweak theory. Say, we take d as incoming and outgoing quark in the left diagram 
of Fig. [U Then, the scattering of W + with the incoming d quark yields a u quark (or c and t if more 
flavors are considered). On the other hand, the crossed diagram on the right in Fig. Q] simply does 
not exist for an incoming d quark, because it is not allowed by the electroweak Standard Model 
couplings. Rather, the incoming quark should be a u quark. In Fig. [3] we display explicitly the 
appropriate pair of Feynman diagrams at leading order. 

Thus, for the partonic forward Compton amplitude (12.171) in the combination t l ^ d = t^ v + t d v 
we effectively sum the contributions of both, crossed and uncrossed, diagrams in Fig. [3] whereas 
for t^y d = fL, — t d v we subtract them. Then, we arrive at the following properties for simultaneous 
transformations /j <-» v and q — > —q, 

t$ d - ^\ (2-25) 
ft* - (-1)^- (2-26) 

Eq. (12.251) implies that the forward Compton amplitude t l ^ d has the same symmetry property as 
in the case of neutral current DIS. For the corresponding coefficient functions and their dependence 
on the Mellin variable n we may repeat exactly the same line of arguments as before leading to the 
conclusion that Ci, Cl (C3) are governed by even (odd) n only. In the other case, Eq. (12.261) shows 
that t l ^ d is antisymmetric under the transformation /j <-> v simultaneously with q — > —q. which 
gives an additional factor (—1) for the l.h.s. of Eq. (12.171) . This alters the Mellin-n dependence 
of the coefficient functions so that we have precisely the opposite assignments, Ci, Cl (C3) being 
entirely odd (even) functions of n only. 

Before moving on, let us briefly comment on the higher order diagrams for charged current DIS 
as illustrated in Fig.|2l For the flavor class fl<i our tree level arguments from above may be literally 
repeated. In the flavor class flqi, on the other hand, crossed diagrams with the same external quark 
flavor do contribute. However, this does not destroy the symmetry properties of the singlet and 
non-singlet combinations. The complete iffi simply sums the crossed and uncrossed diagrams 
and, therefore is still be symmetric under /j <-> v simultaneously with q — > —q, thus Eq. (12.251) 
holds for the u + d combination. In contrast, the contributions from the u — d combination to the 
flavor class /Z02 vanish. The flavor classes fl\\ and //on are excluded in charged current DIS, 
because the flavor changes and, as a consequence, the coupling of one single W + -boson to a quark 
loop is not possible. 

Finally, it remains to relate the coefficient functions C2, C3 and Cl and the OMEs of Eq. (12.151) 
to the Mellin moments of structure functions. To that end, it is convenient to project Eq. (12.31) 
and the analogue of Eq. (12.171) for the hadron forward Compton amplitude T MV onto the respective 
Lorentz structure using projectors 

* s -TFW^"*' (2 ' 27) 
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where all expressions are exact in D = 4 — 2e dimensions. 

With the help of Eqs. (12.271) — (12.291) one arrives at relations between Mellin moments of DIS 
structure functions ( 12.71) , (12.81 ) and the parameters of OPE. On a technical level, this implies a 
Cauchy integration of the analogue of Eq. (12.171) for 7^ v in the complex (O-plane and we recall the 
necessary details in Appendix [Bj 

f dx^- 2 Fi{x,Q 2 ) = Y< C iA n ^^ A L\( n ^ 2 )i ' = 2 > L > ( 2 - 3 °) 
j \ V J 

( dxx n - l F 3 {x,Q 2 ) = Y, C 3j( n ^' a s) A Ll(. n ^) ■ ( 2 - 31 ) 
o i V M J 

To summarize, Eqs. (12.301) and (12.311) provide the basis to obtain Mellin moments of DIS 
structure functions in our approach relying on the OPE and the optical theorem. Furthermore, from 
the careful examination of the symmetry properties of the forward Compton amplitude T^ v and, 
related, the underlying Feynman diagrams, we have deduced corresponding rules for the Mellin 
variable n. In the case of neutral current and charged current vP + vN DIS the structure functions 
F2 and Fl are only functions of even Mellin-« and only functions of odd n for F3. For the case of 
interest in this paper, charged current vP — vN DIS, we encounter only odd functions in n for F2 
and Ft and, only even functions in n for F3, respectively. 



3 Renormalization 

In this Section we briefly recall the necessary steps in renormalizing the operators in the OPE 
and, following the discussion above, we restrict ourselves here entirely to the non-singlet case. 
Starting, say, from the partonic expression (12.171) . i.e. partonic matrix elements of fa, we express 
the renormalized OMEs A J paiton (see Eq. (12.161) ) in terms of matrix elements of bare composite 
operators, 

a, ten =Zm0 o, bare _ (3 1} 

Here and later we suppress other indices and the explicit dependence on n for the operators (12.121) 
The scale dependence of the operator O a is governed by the anomalous dimension y ns , 

^_ «,ren = _^ ^\ (3.2) 

aln^ 

and is connected to the renormalization constant Z ns in Eq. (13.11 ) by 



Ghv 2 -) 



Yns = " T^ Z ns Z ns • (3.3) 
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In order to arrive at explicit expressions for Z ns or Eq. (13.31) . one has to make use of a regular- 
ization procedure and a renormalization scheme. We choose dimensional regularization [24-27] 
in D = 4 — 2e dimensions and the modified minimal subtraction [28,29] scheme, MS. The running 
coupling evolves according to 



a, 



da* 



dln/u 1 An dln/u 2 



(3.4) 



and we have introduced the common short hand notation a s = a s / (4k) . The usual four-dimensional 
expansion coefficients f3 n of the beta function in QCD read p\)= 11 — 2/3 «y etc, with repre- 
senting the number of active quark flavors. The bare and the renormalized coupling, a^ are and a s 
are related by 



(Xj — Z Kj (Xj , 



(3.5) 



where we have put the factor 5 e = exp(£{ln(4ft) — Ye}) = 1 in the MS-scheme and the renormal- 
ization constant Z a reads 



Ax, 



1-— a s + [ rf-^ )«/ + ••• • 



(3.6) 



In this framework, the renormalization factor Z ns in Eq. (13.11) is a series of poles in 1/e, ex- 
pressed in terms of |3 n and the coefficients y^ of the anomalous dimensions from an expansion in 



y(n) = ^ a[ +l y^'^ 

1=0 



( /+i y (/) 

Up to the third order in the coupling constant the expansion of Z ns reads 

Zns 



(3.7) 



1 + as-Yn? + a 2 



+ a: 



{(yS , -23o)(yS ) -Po)yS ) } 



+ ^{3« , -2PoTi', , -2|i 1 T<?} + lTS; 



(3.8) 



The anomalous dimensions yW can thus be read off from the £ 1 terms of the renormalization 
factors at order a s /+1 , while the higher poles in 1/e can serve as checks for the calculation. The 
coefficient functions in Eqs. (12.1 II) . (12.171) . on the other hand, have an expansion in positive powers 
of e, 

oo 

Qns = £ 4 (c g, + za g, + e 2 b g + ...), (3.9) 

1=0 

where i = 2, 3,L and we have again suppressed the dependence on n (and Q 2 /v 2 ). Here Q ns is our 
generic notation for non- singlet contributions obtained for Q a in Eqs. (12.1 II) . (12.171) 
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Due to the presence of 75 in the vertices, the axial-vector coupling in dimensional regulariza- 
tion brings up the need for additional renormalizations to restore the axial Ward-identities. This is 
extensively described in the literature and for the associated renormalizations we use the prescrip- 
tion of Ref. [30, 31] based on relating vector and axial- vector currents. The necessary constant 
Za for the axial renormalization Z5 and the finite renormalization due to the treatment of 75 in the 
MS-scheme are known to three loops [30, 31]. 



The actual calculation of the anomalous dimension (13.21) and the coefficient functions Q ns in 
perturbative QCD proceeds as follows. Using the Lorentz projectors (12.271) — (12.29b we obtain from 
the forward partonic Compton amplitude Eq. (12.171) the partonic invariants 

tip* = If fa, / = 2,3,L, (3.10) 

see also Eqs. (1BI1|) . (|BI2|) . These invariants can be written in terms of the bare operator matrix 
elements as 

ti,ns(x,Q 2 ,a s ,E) = (3.11) 
where i = 2,3, L and the l.h.s. of Eq. (13.1 II) is renormalized by substituting the bare coupling 



constant in terms of the renormalized, see Eq. (13.51) . The wave function renormalization for the 
external quark lines is an overall factor on both sides of the Eq. (13.1 II) and drops out. The terms 
0(p 2 ) on the r.h.s. of Eq. (13.1 II) indicate higher twist contributions, which we neglect. 

Starting with the partonic invariant f ; ns from Eq. (13.1 II) . the renormalization constants Z ns and 
the coefficient functions Q >ns are calculated using the method of projection developed in Ref. [32], 
which consists of applying the following projection operator, 



¥n[f(j>,q)] 



q {n... q Mn} an 

2n~l dpn...dpi*," P,q > 



(3.12) 

P =o 



to both sides of Eq. (13.1 II) . Here q^ 1 ■ ■ ■ is symmetrical and traceless, i.e. the harmonic part of 
the tensor q^ 1 ■ ■ •(f n . 

On the r.h.s. of Eq. (13.1 II) . it is obvious, that the n-th order differentiation in the projection 
operator T n singles out precisely the n-th moment, i.e. the coefficient of 1 fx 11 . All other powers 
of 1 jx vanish either by differentiation or after nullifying the momentum p. The operator fp„ does 
not act on the renormalization constant Z ns and the coefficient functions on the r.h.s. of Eq. (13.1 II) 
as they are only functions of n, a s and £. However, T n does act on the partonic bare OMEs 
^parton' where the nullification of p effectively eliminates all but the tree level diagrams Aj^^, 
because any diagram with loops becomes a massless tadpole and is put to zero in dimensional 
regularization. Finally, the 0(p 2 ) terms in Eq. (13.1 II) . which denote higher twist contributions, 
become proportional to the metric tensor after differentiation. They are removed by the harmonic 
tensor q^ 1 ■ ■ -q^ . On the l.h.s. of Eq. (13.1 II) . T n is applied to the integrands of all Feynman 
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diagrams contributing to the invariants fyparton- The momentum p is nullified before taking the 
limit £ — »■ 0, so that all infrared divergences as p — »■ are dimensionally regularized for individual 
diagrams. This reduces the 4-point diagrams that contribute to t vy! to self-energy type diagrams 
(2-point-functions) accessible to reduction algorithms such as Mincer [33] (see the SectionU). 

To summarize, we find after application of the projection operator T n to Eq. (13.1 II) 

ti,ns(n,^-,a s ,e\ = ^ n ti,ns(x,Q 2 ,a s ,e) (3.13) 

= Q,ns ( n i "^T? C s ,8 J Z ns ( (X s , — J A Q (n,8), 

V r J \ e / 

where z = 2, 3,L. This is our starting point for an iterative determination of the anomalous dimen- 
sions and coefficient functions via the OPE, since the Ci m (Z ns ) are expanded in positive (negative) 

ns tree 

powers of £ while the OME A q ' does factorize after application of the projector <S n . In a series 
expansion in terms of the renormalized coupling a s at the scale fj 2 = Q 2 we can write 

^ns(n) = (O w ) + ^S( n ) + ^S( ,I ) + ^'S( ,I ) + "-) i4 S , ' t,Be(,l) ' (3J4) 

with i = 2, 3, L and recall that we use a s = a s / (4jc) . Then we normalize leading order contribution 
as follows, 

0n) = l, i = 2,3, and ^ 8 (n) = 0, (3.15) 

where the OME A q ' (being a constant) has been absorbed into the normalization of Eq. (I3.15I) . 
With the normalization (13.151) one has 

4%(n) = l, f = 2,3, and cg s (n) = 0. (3.16) 

At first order in a s , expanding up to order £ 2 and suppressing the n-dependence from now on for 
brevity, we find 

'£ = ^ , + c< 1 »+ea« + eX' n ) s + o(e'), 1 = 2,3, (3.17) 

C = 4!L + 4!L + e 2 *i'L + o(e 3 ). (3.18) 
Performing the expansion at a 2 up to order £ we arrive at the equations: 

<S = i{(lS?-Po)lff} + 5{lfl , +2 C g>^} (3.19) 
+ + e + C * } + o (e 2 ) , i = 2, 3 , 

«s = |{4:Lvsv^+4:LYS ) +e{4 2 L+e*}+°( e2 )- 
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Finally for the third order we obtain 

<S = ^{(lS , -2B.)(lffi>-h) 1 ffi>} (3.20) 
{ 3T«°» Tfi» - 2ft, Tfl' - 2fc Tff + 3c« (y£> - ft,) iff } 
+1 {2 T <? + 3c« Y « + 6cf > y<? + 3a<'» (■»£> - ft,) y£> } 

4 Kl + lfl> + a-ffiiff + C (iff " to) lff} + (e) , i = 2, 3 , 

'2. = ^{4i(TS»-to)yS»} (3.2!) 

+s{4:, , , s ^ , + 24 2 , ) , s ^ +a ^:L(^-to)ys , } 

4 HI +4!W»P + +»& (iff - to) iff} + o(e) . 

Eqs. (I3.17I) - (I3.21I) hold for both, even and odd Mellin moments alike and we did not distinguish 
these in our notation. However, from the discussion of the preceding Sections it is clear that the 
respective anomalous dimensions yO and coefficient functions cf^ s describe different physical 

f2) (2) f2l 

processes. In fact, it is well-known that starting from y'- 1 -' and c- ns (and a- ns , £- ns , etc.), they 
differ. The new results of the present paper from Eqs. (13.201) . (13.211) at third order in a s consist of 

(3) (3) (3) 

odd Mellin moments for c 2 4 and c L m and even moments for c 3 ^. Below in Sec. @] we present 
numerical results for them, while complete expressions through rational numbers are deferred to 
Appendix [Al 



4 Calculation 

In the previous Sections, we have laid the foundations to our calculation of Mellin moments of the 
DIS charged current structure functions F 2 p ~ vP , F^ p ~ vP and f^ p ^ vP together with their respective 
coefficient functions and anomalous dimensions. To that end, following Refs. [1,2, 8,9], we have 
calculated the Lorentz invariants of the parton Compton amplitude tf^, I = 0, 1 , 2, 3, i = 2, 3, L, as 
given in the l.h.s. of Eqs. (I3.17I) - (I3.21I) from first principles. All contributing Feynman diagram 
were generated and then projected by one of the Lorentz projection (12.271) — (12.291 ). Subsequently, 
the application of Eq. (13.121) for the harmonic projection <S n extracts all contributions to the given 
Mellin moment, which are finally solved in terms of rational numbers, values of the Riemann zeta 
functions and SU(N C ) color coefficients Ca, Cp and tif. Due to the large number of diagrams 
involved in the calculations up to order oc^ sufficient automatization is necessary. Therefore the 
calculations are organized in detail as follows: 

• All Feynman diagrams are generated automatically with the program QGRAF [34]. This pro- 
gram generates all possible Feynman diagrams (and topologies) for a given process in some 
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Lorentz 
invariant 


Structure 
function 


tree 
o(a°) 


one-loop 

o(«J) 


two-loop 

o(aj) 


three-loop 

o(aj) 


sum 




r VP-VP 


1 


4 


55 


1016 


1076 


^L,ns 


r vP-vP 
L 


1 


4 


55 


1016 


1076 


^3,ns 


^3 


1 


4 


63 


1246 


1314 


in total 












3466 



Table 1 : The number of diagrams involved in the calculation of the vP — vP charged current DIS 
structure functions F<i, Fi and Ft, at tree level, one-loop, two-loop and three-loop, respectively. 

special format. The program works very effectively, producing a database with thousands 
of diagrams within seconds only. For charged current DIS we have obtained from Qgraf 
2, 10, 153 and 3468 diagrams for the tree, one-loop, two-loop and three-loop contributions, 
respectively. 

• For all further calculations we have relied on the latest version of the symbolic manipulation 
program Form [35,36]. For the further treatment of Qgraf output, such as analysis of the 
topologies, the explicit implementation of Feynman rules etc. we have adapted a dedicated 
Form procedure conv.prc from previous work, e.g. Ref. [9]. Most importantly, this proce- 
dure tries to exploit as many symmetry properties of the original Feynman diagrams in order 
to reduce their total number. The upshot of these efforts are presented in Table[T]order by or- 
der for structure function corresponding to the different Lorentz projections. As one can see, 
the number of diagrams obtained for F^ p ~^ p is always larger than for F^ or F^ p ~^ p . 
The reason is that in the former case we can not apply certain symmetry transformations due 
to the presence of 75 in the vertices. The database for i 7 ^ p ~^ p and F^ p ~ produced by us 
does almost coincide with the one used in Ref. [16], except for small modifications. The 
database for F^ p ~ is completely new. 

• For the calculation of the color factors for each Feynman diagram we have used the Form 
package color.h [37]. 

• The actual calculation of the Mellin moments of the Feynman integrals has made use of 
Mincer. The detailed description of this program can be found in Ref. [33] for the Form 
package mincer.h. For organization of the work of (a slightly modified version of) Mincer 
with the input databases we have used a dedicated database program MlNOS [14, 15]. 

• Finally, on top of Mincer and Minos some shell scripts managed the automatic runs of 
both programs for different parts of the calculation. This facilitates the bookkeeping of 
different input parameters for MlNOS and Mincer due to different Lorentz projections, 
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orders of a s etc in distributed running. Moreover, the shell scripts also organized the final 
summations over the flavor classes as well as the output of all final results. 

Next, let us discuss the various checks we performed on the results of our calculations. First 
of all, we have tested our set-up by a recalculation of some known even Mellin moments for F% 
Fl and odd moments for F3 to find agreement with the published results of Refs. [9, 14-16]. Then, 
most importantly, we have checked gauge invariance, i.e. we calculated all our diagrams for all 
results for all Mellin moments presented in this article with a gauge parameter £, for the gluon 
propagator, 

T 

We kept all powers of £, (up to ^ 4 in three loops for this calculation) through the entire calculation. 
Since parton structure functions are physical quantities any dependence on the gauge parameter ^ 
must disappear in the final result. This was indeed the case after summing all diagrams in a given 
flavor class. Furthermore, we have compared the anomalous dimensions y ns Eq. (13.21 ) as calculated 
by us from Eqs. (13.17lHl3.21l ) up to three loops with the results available in the literature [1] and 
found complete agreement for both, even and odd Mellin moments. In addition, the coefficient 
functions for the structure functions pv p - vP ^ { = 2, 3,L at two loops are known from earlier work. 
Our two-loop results as obtained from Eq. d3.17lHl3.20l) coincide with Refs. [3-7]. 

Finally, let us mention a few words on the hardware requirements. All calculations are CPU- 
time and disk space consuming, especially for the higher Mellin moments (higher n values). They 
were typically performed on an 64-bit AMD Opteron 2.2 GHz Linux machine with 4 GByte of 
memory. For example, the calculation of ?3 ins for n = 10 took 56 days with the gauge parameter 
included, while the calculation of both, ?2,ns and fL.ns for n = 9 on the same machine needed 33 
days. For comparison, the calculation of lowest Mellin moment n = 1 for both projections ?2,ns and 
fL.ns consumes less than a hour, whereas for n = 2 for ?3. ns one needs approximately a couple of 
hours, always with the full gauge parameter dependence. At intermediate stages the calculations 
required also a large amount of disk space. Although the programs calculate diagrams one by 
one at the time the size of intermediate algebraic expressions for some diagrams can grow up to 
20 GByte of a disk space (for instance for n = 10 three-loop diagrams). On the other hand, the 
final result for any of the Lorentz invariants occupies some KBytes only. With access to improved 
hardware, we plan to push the calculation further to n = 16, cf. Ref. [18]. 



5 Results 

Following the steps outlined above in Sec. |4] we arrive at the results for the coefficient functions 
C2,ns, ^-L,ns at the odd-integer values n = 1, . . . , 9, and for C3 ins at the even-integer values n = 
2, . . . , 10, up to order oc^. The third order expressions represent new results of this article. Using 
a s = a s / (4k) and the shorthand notation for the n-th moment Cf* n = Ci jns (n) we find the following 
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Cf 3 = 1.333333333 a s + a 2 (52.40384466- 3.925925925 n f ) (5.7) 



numerical values at the scale /j r = /jf = Q, 

C 2 ns ! = 1, (5.1) 
Cf s 3 = l+3.222222222a s + a 2 (72.32720288- 11.125/2/ (5.2) 

+ a 3 s (1948.031519-496.5427343/^ + 14.20173594 /z/) , 
Cf 5 = l+8.725925925a s + a 2 (220.4151827 -22.64048559/i/ (5.3) 

+ a 3 & (6925.814438- 1347.125829^+32.94421923 n/), 
Cf p = 1 + 13.43677248 a s + a 2 (386.291 1104- 33. 10212484/i/ (5.4) 

+ a] ( 1 3505 . 1 6600 - 2298 . 472900 n f + 52 . 34745652 n f 2 ) , 
Cf 3 = 1 + 17.47820105 a s + a 2 (555.27201 17 -42.50367619/i/ (5.5) 

+ a] (20990.73668 - 3278.689323 n f + 71 .3 1040423 n f 2 ) , 
Cf A = +2.666666666a s + a 2 (61. 33333333 -4.740740740 /i/ (5.6) 
+ a 3 s (2313.911655 -405.2001359/1^ + 10.20576131 n/) , 

n f ) 

+ a] (2584. 178446 - 406.0509532 n f + 1 1 .5967078 1 n f 2 ) , 
C™ 5 = +0.8888888888a s + a 2 (44.23466187-3.012345679^) (5.8) 

+ a* (2451.068575 - 360.6487058/^ + 10. 15089163 « f 2 ) , 
Cf p = 0.6666666666a s + a 2 (38.25090234 - 2.440740740 /i/ (5.9) 

+ a 3 & (2290.679208 - 321.7773285/^ + 8.868930041 n/) , 
Cf 9 = 0.5333333333a s + a 2 (33.82305394 -2.05671 9576 /i/ (5.10) 

+ «3 (2146.302724 - 290.9906309/^ + 7.868039976/1/) , 
Cf 2 = 1 -1. 777777777 a s + a 2 (-47.07704646 -0.09876543209/?/ (5.11) 

+ ^(-2359.001407 + 305.2538856/^-6.864103442/1/), 
Cf A = l+4.866666666a s + a 2 (90.15322509- 1 3.25902469 /i/ (5.12) 

+ a^ (1478.747872 -491.0449098^+ 11. 77903924 nf) , 
Cf fi = l + 10.35132275a s + a 2 (258. 8595696-25. 14210054/i/ (5.13) 

+ a^ (7586.717646 -1458.855783^+ 32.73965909 n/), 
Cf s 8 = l + 14.90026455a s + a 2 (433.2106396-35.58166191/i/ (5.14) 

+ a 3 s (14862.60949 -2469.591886/1^ + 53.25812942/1/) , 
Cf s 10 = 1 + 18.79152477 a s + a 2 (605.9424494 -44.87506803 /i/ (5.15) 

+ a 3 & (22806.38215-3482.933316^ + 72.83344233/1/). 

Exact analytical expressions of these moments also with complete dependence on the color coeffi- 
cients are given in Appendix lAl 

As was mentioned before, the two-loop coefficient functions in Eqs. (I5.1I) - (I5.15I) agree with 
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the results in Refs. [3-7]. In addition, Eq. (15.11) for is nothing else but a manifestation of the 
Adler sum rule for DIS structure functions, 

l 

/ f (^fr^-^frfi 2 )) =2, (5.16) 
o 

which measures the isospin of the nucleon in the quark-parton model and does not receive any per- 
turbative or non-perturbative corrections in QCD, see e.g. the discussion in Ref. [38]. Therefore, 
Eq. (15.11) is another important check of the correctness of our results. 
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Figure 4: The first ten integer Mellin moments of the third order non-singlet coefficient functions 

(3) (3) 

c\ „ s (left) and c L ni (right) for charged current DIS with rif = 4 flavors. For the even moments, the 
flavor class //02 has been omitted, i.e. //02 = (see also the discussion in Sec. [2]). 



For illustration, let us plot the Mellin moments of the coefficient functions c\ c\ L and c, L 
at three loops. The new non-singlet results (blue diamonds in Figs. [4] and [5]) exhibit a similar 
smooth pattern as the known results (red squares in Figs. |4] and [5]). Thus, it is feasible to use these 

(3) 

moments for an approximate analytic reconstruction of the yet unknown coefficient functions c\ ns 
for (yP — v.P)-DIS prior to "all-n" calculation and similar to e.g. Refs. [36,39]. We will do so in a 
companion paper [21]. 

Furthermore, we see that the respective values for odd and even moments, for instance on 
the left in Fig. |4] do confirm that differences between C2^ P+ ^ P and c^l^ P ^ P are numerically 
small. This observation (see Fig. |5]) provides also a posteriori justification for the extrapolation 
procedure from odd to even moments for C3 in Refs. [40,41]. There, available information on odd 
moments [16] used in fits of CCFR data [42] to the structure function xFj, at NNLO in QCD and 
beyond. A further discussion of this and related issues is given in Ref. [21]. 
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Figure 5: The first ten integer Mellin moments of the third order non-singlet coefficient functions 

(3) 

c 3 ns f° r cnar g e d current DIS with rif = 4 flavors. For the odd moments, the flavor class /Z02 has 
been omitted, i.e. /Z02 = (see also the discussion in Sec. El). 



6 Conclusions 

In the present paper we have presented new results for Mellin moments of the charged current 
DIS structure functions P-uP , F^ p ~ yP and F^ p ~^ p including the perturbative QCD corrections 
to three loops. In the former case (/% Fi) we have computed the first five odd-integer Mellin 
moments while in the latter case (F$), the first five even-integer moments have been given. Our 
efforts are part of an ongoing program [1-9, 14-18] to calculate perturbative QCD corrections of 
all DIS structure functions to three-loop accuracy. 

Within the framework of the OPE and the optical theorem we have calculated Feynman dia- 
grams of the forward Compton amplitude 7^ v in Mellin space. In our presentation we have em- 
phasized the symmetry properties of and their relation to charged current vP±\P DIS, which 
was a crucial point in setting up the databases of Feynman diagrams. We have performed various 
checks on our computation. Most prominently, we have kept all powers of the gauge parameter 
£, throughout the entire calculation to check that any ^-dependence vanishes in our final results. 
Furthermore, we agree with the literature as far as the two-loop coefficient functions [3-7] and the 
three-loop anomalous dimensions [1] are concerned. 

The discussion of phenomenological consequences of our Mellin space results is deferred to 
Ref. [21]. Future research will be devoted to the calculation of some higher Mellin moments, 
potentially n = 11, . . . , 16, depending on the available hardware infrastructure. Subsequently, we 
will also focus on an "all-n" calculation in Mellin-n space with methods of Refs. [1,2, 8, 9], since 
all databases for Feynman diagrams contributing to F^ p ^ p , F^ p ^ vP and f^p-vp ^ ava {i a bi e now 

Form files of these results can be obtained from the preprint server |http: //arXiv. org] by 
downloading the source of this article. Furthermore they are available from the authors upon 
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request. 
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A Appendix 

In this Appendix we present the analytic expressions up to order for the coefficient functions 
C™, C" s at the odd-integer values n = 1, . . . , 9 and for C" s at the even-integer values n = 2, . . . , 10. 
The notation follows Sec. [5] with Ca and C F being the standard QCD colour factors, Ca = N c = 3 
and Cf = (Af? — 1)/ (2iV c ) =4/3, and stands for Riemann's ^-function. The coefficient functions 
for the structure function F2 at the scale ji r = jif = Q are given by 

cji = 1, 61i) 

Cf 3 = 1 CS2) 



^ 29 
+ a s C F — 

+ alC F nj-(^- 

+ alC F 2 {- 



267 

~1$2~ 
217235 

10368" + 36 ^ 3 



/25855 



+<QC F (^32--43C3 

2 /641563 100,, \ 

+ <C F «/( 1 C3 

s F / V 69984 81 s ) 

3 2 / 35314337 1577., 50. 

+ a * CF n H--699840""^ 3 + T^ 4 

3 3 /57093841 103235 55 1520. 

+ asC ^ I 373248 + ^24~^ + T^ 4 "^ 5 

3 / 132494393 39203,, 50,, \ 

+ a a C A C F n f (- 69984Q + — ^3 - y^J 

3^2/ 614328541 16171 „ 55 „ _ 

+ a C 4 Cp 1 C3 C4 + 40^5 

s A F \ 2799360 45 s 2 s w 



3 2 / 490358569 344929 55 „ 1070. 

+ ° sQ C H~699840 405~^ 3 + ^ 4 + ^-^ 

Cf 5 = 1 ©3) 



„ 589 
+ a 2 C F n f (^- 



2750819 
162000 
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. 2 / 30297101 cny \ 

2 / 35848409 312 
+ a s L A t F { 324()00 — c, ; 

3 2 / 10958051 728 
+ a * C F n f { 486000 + 405 



3 ^ 2 /10958051 728, \ 

+ a ^ n f 1^86000- + 405^ 

, 2 / 38720716199 257668. 364 \ 
+ "A" 255150000 -^725"^ 3 + l5" C4 J 

3 3 /132383443837 5617861 1414. c , OJ , \ 

+ I 2187000000 + ^0l2T C3 + ^ 4 ~ 568C5 J 

3 / 50351664421 2187697. 364. \ 

+ a s C A C F n f [- 122472000 + "^4175-53 - -^J 

+ alC A C/ ( 9620298263 + ^^3 - ™ - 108&) 
s A F V 63787500 23625 S3 25 S4 S V 
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+ aiC A l C F [ C3 H C4 + 56OC5 , 
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, / 94605 2961 283 _ 1810 

+<*sC A c F [ 4898880000 2T^3 

, 2 /3438632355495191 4180,, 
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s F f \ 68052791520000 1701 b 
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The coefficient functions for the structure function Ft at the scale /j r = /jf = Q are given by 

Cl\ = a s C F 2 dAl6) 
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B Appendix 



In this Appendix, we recall a few technical steps necessary to arrive at the relations (12.301) . (12.311) 
between the parameters of OPE and the Mellin moments of DIS structure functions. To that end, 
we would like to put particular emphasis on the symmetry properties of the hadron forward Comp- 
ton amplitude T^(p,q) in Eq. (12.101) under the transformations /j <-> v and q — > —q. Since our 
discussion in Sec. [2] was largely based on Feynman diagram considerations at parton level it re- 
mains to link the line of arguments to the OPE of Eq. (12.1 II) . to the analogue of the OPE (12.171) for 
7^ v and, eventually to the Mellin moments of Fz, Ft, and Fi in Eqs. (12.71) . (12.81) . 

Let us start by observing, that the OPE (12.1 II) for 7^ v gives rise to a series expansion in terms 
of co similar to Eq. (12.171) . This series being valid for unphysical CO = \/x — > only (recall the 
Bjorken variable < x < 1) is related to the physical Mellin moments of F2, F3 and Fl by means 
of a Cauchy integration. Here, the behavior of 7^ v under the mapping co — > -co (q — > —q) becomes 
relevant. Applying the Lorentz projectors (12.27b — (12.29b to we obtain 

Ti((0, Q 2 ) = /f i>(G), Q 2 ) = 2 J> n C LJ L, %,cj) A J nucl {n^ 2 ) , i = 2,L. (El) 

n,j \ V J 

In Sec. [2] we discussed which coefficients Clj{u) survive in the OPE (12.1 II) . As consequence 
the sum (IB! II) runs over even n for the neutral current structure functions F2, Fl and the charged 
current (singlet) structure functions F~ + , F^ + and therefore 7](— CO, Q 2 ) = Tj((0,Q 2 ). For 
the charged current (non- singlet) structure functions F~ ~ and F£ ~ we sum over all odd n. 
Therefore 7}(— CO, Q 2 ) = — 7/(co, Q 2 ). Furthermore, we have 

T 3 ((o, Q 2 ) ee if rpv(co, 2 2 ) = 2j>"C 3j (w, ^,a,") A^ ucl (n,^ 2 ) , ©2) 

where we sum over all odd n for the neutral current F3 and charged current f^ p+vN . As conse- 
quence we have obviously T3 (— (0, <2 2 ) = —73(0), Q 2 ). On the other hand, for the charged current 
pvP-vN sum mns oyer eyen ^ an( j f^gj-gfoj-g 7 3 (_(j) 7 g 2 ) = 73(00, <2 2 ). 

In applying the Cauchy integration to both sides of Eqs. (IB 111) . (IBI2I) with a contour C as shown 
in the Fig. [6] we exploit the fact that the Lorentz projected forward Compton amplitude is an 
analytic function of complex variable CO. The branch cuts extend along the real axis for co < — 1 
and co > 1 because of kinematical constraints from Bjorken x and symmetry properties. We divide 
both sides of Eqs. (IBI1I) . (IBI2I) by 2ni(D m and pick up the appropriate residues on the r.h.s according 
to 

1 f co" 

— <b d(0— = 8 n , ra _i . <E3) 
2m / co"' 



c 



For the l.h.s. one obtains 



1 ie +©o+i£ 1+ie — °o— ie 



2*j^ t ^ q » = mI + y + y + J + y 



C +oo-ie +l+ie -°o+i£ — 1— ie 
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Figure 6: The contour C of the Cauchy integration in the complex (O-plane. 
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J —JmTi((0,Q 2 ), i = 2,3,L, 



l 



where we have used that for the physical cases with either even or odd n in Eqs. (IBI1I) . (IBI2I) the 
whole combination (rfro/co'") 7/(co, <2 2 ) occurs to be symmetric under (0 — > -co, because of the 
symmetry properties of 7} and the restriction for m Eq. (IBI3I) . Thus, 

^ 7}(ffl, 2 2 ) ^ +^ 7i(ffl, 2 2 ) , » = 2, 3,L, (E5) 

and, in addition, we have used taken the discontinuity across the branch cut 

Im^((0,e 2 ) = ^(7; ((o + ie,<2 2 )-7;((o-ie,e 2 )) , i = 2,3,L. dEJ6) 

Changing variables (0 — > 1/jc in Eq. (IB 1 41) and using Eq. (IBI3I) for the r.h.s. of Eqs. (IBI1I) . (IB 121) one 

gets 

i / dxx"- 1 hnTi(x, Q 2 ) = Y,Qj (n, ^-,a s ) A{ ud (n, a/ 2 ) , / = 2, 3,L. ©7) 

Q J \ M J 

Applying the Lorentz projectors (I2.27l) - (|2.29b to the optical theorem Eq. (12.91) and using Eq. (12.31) 
for hadronic tensor we obtain 

-ImT^Q 2 ) = -Fi(x,Q 2 ), i = L,2, dB8) 

7t 

-Imr 3 (x,2 2 ) = F 3 (x,<2 2 ). (E9) 

The combination of Eq. (IB[7I) with Eq. (IB[8I) concludes the derivation of Eqs. (12.301) . (12.311) with 
special emphasis on the transformation properties of under 0) — > — 0) (<? — > — 
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